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Abstract. We use the perturbative renormalization group to study classical 
stochastic processes with memory. We focus on the generalized Langevin dynamics 
of the (f)'^ Ginzburg-Landau model with additive noise, the correlations of which are 
local in space but decay as a power-law with exponent a in time. These correlations 
are assumed to be due to the coupling to an equilibrium thermal bath. We study both 
the equilibrium dynamics at the critical point and quenches towards it, deriving the 
corresponding scaling forms and the associated equilibrium and non- equilibrium critical 
exponents 77, i^, z and 9. We show that, while the first two retain their equilibrium 
values independently of a, the non-Markovian character of the dynamics affects z and 6 
for a < ac{D,N) where D is the spatial dimensionality, N the number of components 
of the order parameter, and ac{x, y) a function which we determine at second order in 
4 — D. We analyze the dependence of the asymptotic fluctuation-dissipation ratio on 
various parameters, including a. We discuss the implications of our results for several 
physical situations. 



Critical dynamics with correlated noise 



2 



1. Introduction 

For more than 30 years critical dynamics have been explored with field theoretical 
methods [Il|2l[3llll[5l[6l[71[8l|9]. A variety of dynamic models were introduced to describe 
the collective evolution of systems close to critical points. Among these, the most 
common ones are the dynamics of non-conserved or conserved order parameters which 
successfully describe the evolution of uniaxial magnetic systems close to the Curie point 
or the dynamics of binary alloys close to the demixing transition, respectively. These 
problems as well as many of their generalizations discussed in [H |5l [71 E] are classical in 
the sense that their stochastic nature can be essentially ascribed to thermal fiuctuations. 
Therefore, as in the simpler diffusion processes, the evolution of the interacting degrees 
of freedom — described via a field — is modeled by a functional Langevin equation 
in which the coupling to the environment is responsible for both thermal fiuctuations, 
encoded in a stochastic external noise, and the friction. If the environment, which acts as 
a thermal bath for 0, is in equilibrium at a certain temperature the space and time 
dependence of the friction coefficient and of the correlations of the thermal fiuctuations 
are related via the fiuctuation-dissipation theorem. The selected model bath determines 
then the remaining functional form of the noise-noise correlation and the usual choice 
is to take it to be delta-correlated in time which corresponds to white noise. 

In cases in which the initial conditions of the system are drawn from an equilibrium 
Gibbs-Boltzmann distribution, or the system is allowed to evolve for a sufficiently long 
time such that this distribution function is reached, the space-time behavior of dynamic 
quantities is characterized by scaling laws in which the usual static exponents (z/, 77 
etc.) appear but a new critical exponent z is needed to relate the space and time 
dependencies. Besides the analysis of equilibrium dynamics, field-theoretical methods 
allow one to study the nan- equilibrium dynamics after a sudden quench from a suitably 
chosen initial condition to the critical point ^U\. A Gaussian distribution of the initial 
field configuration with zero average and short-range correlations mimics a quench 
form the disordered phase [TOl. A distribution with non-zero average but still short- 
range correlations describes (in the case of a scalar field 0) a quench from the ordered 
state [m [12] . In these non-equilibrium cases a new critical exponent — usually denoted 
by 6 and called the "initial slip exponent" — characterizes the short-time behavior of 
the average order parameter as well as of the correlation and response functions [10] 
(see, e.g., [I3l [HI [I5] for summaries). 

In all the studies mentioned above, the noise is assumed to have a Gaussian 
distribution (as a consequence of the central limit theorem) with no temporal 
correlations, i.e., to be Gaussian and white. However, the coupling to thermal reservoirs 
yields, in general, non-Markovian Langevin equations in which the noise is correlated 
in time and the friction coefficient has some memory [T6l \T7\ [T8] . In many situations 
of practical interest a full treatment of such a colored noise is therefore necessary. For 
example, the escape rate of particles confined within a potential well crucially depends on 
the statistics of the thermal bath [T9t [20l [2T] as observed in the desorption of molecules 
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from a substrate undergoing a second-order phase transition which effectively provides 
a colored noise for the stochastic dynamics of the molecules (this phenomenon is called 
Hedvall effect and it was studied theoretically in, e.g., [22]). Another important instance 
is the stochastic Burgers modeling of turbulence where the noise is correlated in both 
time and space [23]. This equation is, in addition, closely related to the Kardar-Parisi- 
Zhang description of surface growth that was analyzed with spatially correlated noise 
in, e.g., [2H |25] and references therein. The so-called fractional Brownian motion [26] 
is actually defined by the non-Markovian nature of the process, which can be traced 
back to the temporal correlations characterizing the thermal noise. The physical 
circumstances in which a temporally correlated noise arises are manifold including 
polymer translocation through a nanopore [27| or the effective description of a tracer in 
a glassy medium [28]. A review of the effects of colored noise in dynamical systems is 
given in [29]. Last but not least, the environment fluctuations in quantum dissipative 
systems give rise rather naturally to temporally correlated contributions |18] . 

In the present study we explore the influence of colored noise on the equilibrium 
and non-equilibrium critical dynamics of the A^-component 0^ Ginzburg-Landau model 
with 0{N) symmetry and non-conserved order parameter. We focus on the case of noise 
correlations that decay as a power-law in the time-lag, such that no time scale can be 
associated to it at least for sufficiently long times. To our knowledge critical dynamics 
with such a colored noise have not been theoretically investigated yet. Our first aim is 
to present a general equilibrium analysis that allows us to determine the classes of noise 
which affect the critical relaxation. We use the perturbative renormalization group (RG) 
technique to calculate the critical exponents, which turn out to be modified for a certain 
type of noise correlations only. Our second aim is to understand the non-equilibrium 
dynamics in the presence of such kind of colored noise. A short time after the critical 
quench, the dynamics of the system can be described in terms of the non-equilibrium 
critical exponent 6 which depends on the properties of the colored noise. In the critical 
coarsening regime we define an effective temperature [30l [31] from the long-time limit 
of the so-called fluctuation-dissipation ratio [c.f., Eq. f l4.40p for its definition] which was 
proposed to be an additional universal property of the non-equilibrium dynamics |32j . 

was calculated for a variety of critical processes with white noise [HI [33l |34] • We 
show that turns out to be affected by the color of the noise and we compare our 
results to the corresponding ones for the fractional Brownian motion [28] . 

The paper is organized as follows. In Sec. 2 we present the model and we obtain 
the associated dynamic action by using the path-integral formalism [H IH |2l [TJ El |35] . 
We define the response and the correlation functions, we present their scaling behavior 
and we define the critical exponents. In Sec. 3 we focus on the equilibrium properties 
of the model, we calculate the propagators in the frequency and momentum domain, 
presenting the Feynman rules of the perturbation theory. This section contains details 
on the perturbative RG analysis which allows us, on the one hand, to determine the 
type of noise which alters the critical dynamics and, on the other, to calculate the 
dynamical critical exponent z up to second order in e = — D where is the 
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upper critical dimensionality, which we find to be 4. In Sec. 4 we focus on the non- 
equilibrium dynamics. We derive the critical scaling behavior of the linear response 
and correlation functions and we calculate the non-equilibrium critical exponent 6 (the 
"initial slip exponent", see [10]) up to first order in 4 — D. Using the scaling forms of 
the correlation and linear response functions we deduce the effective temperature after 
a critical quench. Finally, in Sec. 5 we summarize our results and we discuss several of 
their possible applications together with plans for further investigation. The details of 
the calculations are reported in the Appendices. 



2. The model 



In D spatial dimensions we consider an iV-component non-conserved order parameter 
4>{x) whose static behavior near the critical point is governed by a Hamiltonian of the 
Ginzburg-Landau type 



n\(j)]= I d^x 



i(V0T + lr0^ + |, 



(2.1) 



g is the strength of the non-linearity that drives the phase transition, r is the control 
parameter for it, and the coefficient in front of the elastic term oc (V0)^ has been 
absorbed in the definition of the field. The time evolution under purely dissipative 
dynamics (model A according to the standard classification of [^) in the presence of 
memory is described by the Langevin equation 

dt' T{t - t')dt'(l){x, t') + ^ ^ = C(x, t), (2.2) 

-T 6(f){x, t) 

where —T is the initial time of the process and C is a zero-mean Gaussian colored noise 
with 

{ax,t)Q{x',t')) = r'r{t -t')6{x - x')6,,. (2.3) 

r is a positive and symmetric kernel, i.e., r(t — t') = r(t' — t) > 0. Depending on the 
specific system of interest, in the minimal dynamical model in Eq. (12. 2p one might need 
to include an inertial term Md^cj). An order of magnitude analysis suggests that such 
a term is relevant for times shorter than M/ J dt T(t), whereas the effects of inertia 
can be neglected for longer times, indeed the regime we are interested in. For this 
reason we shall omit this term from the outset as long as the contribution of friction 
does not vanish. Note that the function F determines both the noise-noise correlation 
[Eq. (12. 3p ] and the time-dependent retarded friction coefficient [Eq. (12. 2p ] since we have 
assumed the thermal bath (which is weakly coupled to the system) to be in equilibrium 
at temperature [we set ks = 1]. The Markovian examples of this dynamics are 
characterized by a (5-correlated ('white') noise, i.e., 

F(t) = 27^5(t), (2.4) 

where 7^ is the friction coefficient. In this case Eq. (12. 2p has been extensively studied 
both in and out of equilibrium, see, e.g., [31 HDl [HI [36] • Such Ohmic dissipation is the 
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simplest form of short-range correlated noise. It can be formally obtained as the limit 
to — of the exponentially correlated Ornstein-Uhlenbeck (OU) process 

rou(t) = 2^e-l*l/\ (2.5) 
to 

where the finite characteristic relaxation time to plays the role of an internal scale. Under 
renormalization one expects the exponentially correlated noise to become equivalent to 
a white (delta-correlated) one, Eq. (12. 4p . and the critical behavior of the OU process be 
identical to the Markovian one. In the absence of an internal scale, instead, there is no 
reason to expect a white noise limit and the critical behavior might be affected. The 
simplest example with no explicit time scale is 

T(t) = — --^ with a > 0. (2.6) 

Te[1 - a) 

We denote by F^; the Euler F-function to avoid confusion with the noise kernel F. For 
a > 1, i.e., super Ohmic dissipation, expression (12.61) is not integrable, unless a short- 
time cut-off and thus an internal scale is introduced. One can show that under naive 
scaling (introduced in Sec. 3) the Fourier or Laplace transform of F(t) generate a white 
noise vertex that dominates over the colored noise part. Hence, the appearance of a 
cut-off scale suggests the non relevance of the colored noise for super- Ohmic dissipation, 
i.e., a > 1. This statement will be made precise in the following. Instead, for sub- 
Ohmic dissipation, i.e., a < 1, the noise is truly long-range correlated and its influence 
on the dynamics will turn out to be non-trivial. The naive cross-over value between 
these two cases is a = Oc = 1, that is white noise or Ohmic dissipation. In the presence 
of interactions we shall show that this scenario is slightly modified, with the cross-over 
value ac{D, N) depending upon D and A^. 

A functional- integral representation of the stochastic process, Markovian or not, is 
better suited for an analytic treatment of critical dynamics than the Langevin equation 
(12. 2p . In particular, it allows one to express the average {■■■)(; over the possible 
realizations of the noise ( in Eq. (12. 2 p as a functional integral (which will be denoted 
by (■ ■ ■) in what follows) 

(■ ■ ■)c = I [d0d0] ■ ■ ■ e-^[<^'^l (2.7) 

over and an auxiliary field with S = Sq + Sint — In Vjc [D IS IH [35j|l . 

dt' V{t - t')dt'M^, t') + (r - V2)0.(x, t) 



So = / d^x / dt 0i(f,t) 



T 



and 



//•OO Pt 
d^x J ^tj dt'^i{x,t)T{t-t')'^i{x,t') (2.8) 



Sint= I d^x I dt ^(f)i{x,t)(j)i{x,t)(l)j{x,t)(f)j{x,t). (2.9) 



I In the presence of colored noise no discretization problems arise, see, e.g., [35] . 
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We used Einstein's convention of summation over repeated indices. The zero-source 
functional integral is identical to 1 due to the normalization of the noise probability 
distribution. 'P/c[0(3;, —T)] is the statistical weight of the initial condition. The 
auxiliary field ^ is conjugated to an external perturbation h, in such a way that if 
^[0, h] = ^[(j)] — (j) ■ h, the linear response of the order parameter to the field h is given 
by 



R{x-x';t,t')S, 



6hj{x',t' 



= (2.10) 

h=6 

where {■ ■ is the average over the stochastic process in the presence of the external 
perturbation, i.e., Eq. (12. 2p with "H H- 'H[0, /i]. The response function is causal 
irrespectively of the noise statistics and the Jacobian of the transformation of variables 
from C to which allows us to write the average over the stochastic process as in Eq. (12.71) 
is also a factor with no consequences [35]. In addition to the (linear) response function, 
we shall consider below the correlation function of the order parameter, defined by 

C{x - x', t, t%j = {4>i{x, t)<p,{x\ t')) (2.11) 

where we assumed translational invariance in space. The action Sq + Sint is the sum 
of two contributions each one made of several terms. The part with density (j)i6T-L/6(j)i 
represents the deterministic dynamics whereas the remaining part is due to the coupling 
to the bath. The latter consists of the friction term and the noise-noise correlation and 
both involve the kernel F. In this formulation the problem is recast in the form of a field 
theory in D + 1 dimensions with two vector fields, the analysis of which can be done via 
standard field-theoretical tools, such as the renormalization group (RG) approach that 
we shall use below. 

Since, in general, there is no tractable Fokker-Planck equation for the non-Markov 
stochastic processes we are presently interested in, the usual and relatively simple proof 
of equilibration explained in, e.g., [STJ [38] for the white-noise problem does not apply. 
However, we recall here that Eq. (12. 2p is an effective description of the dynamics of a 
classical system with Hamiltonian "H' which is weakly and linearly coupled to a (large) 
equilibrium bath of harmonic oscillators at temperature acting as a source of the 
stochastic noise ( effectively induced by such a coupling. Indeed, the temperature that 
characterizes the correlations of the noise in Eq. (12. 3p is /3~^, whereas the distribution 
of the frequencies of the harmonic oscillators within the bath determines the functional 
form of r. In addition, F appears in Eq. (12. 2p and Eq. (12. 3p in such a way to ensure the 
fluctuation-dissipation theorem for the bath variables [see, c.f., Eq. (I2.12p ]. As a result, 
even with this effective non-Markov dynamics the system should still lose memory of 
its initial condition and equilibrate with the equilibrium bath of oscillators, resulting in 
a canonical distribution e'^'^^^y Z{(3) of one-time quantities at sufficiently long times 
(possibly divergent with the system size) where is the partition function and T-L 

differs from "H' by a term which is quadratic in the relevant degrees of freedom (see, 

§ is purely imaginary and it is sometimes written as in the literature. 
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e.g., [IHIEH] for details). The asymptotic critical equilibrium dynamics is expected to be 
described by the limit T — oo of the action in which one neglects the specific distribution 
Vic of the initial conditions that in any case should be forgotten dynamically. Since 
we shall be interested in the critical dynamics, we set (3 = 13^ and we absorb this 
constant into a redefinition of the fields and of the coupling constant g. In equilibrium 
the response and the correlation functions are invariant under time translations, i.e., 
R{x, t, t') = R{x, t - t') [see Eq. fl2:T0|) ] and t, t') = C{x, t - t') [Eq. fim]) ] , 
and they are related to each other by the fiuctuation-dissipation theorem (FDT) that 
reads R{x,t) = —f3dtC{x,t)Q(t), where t represents the time delay, 9(t < 0) = and 
G(t > 0) = 1 [Jj], and which is completely independent of the specific characteristics 
of the system and the bath apart from its temperature. (A proof of this relation for 
generic colored noise Langevin dynamics can be found in |35].) Once the latter has been 
absorbed in the redefinition of (pi and g the FDT becomes 

R{x, t) = -dtC{x, t)e(t), (2.12) 

and this is the form that we shall use in our calculations. Moreover, the time-dependent 
correlation is invariant under time- reversal, i.e., C{x,t) = C{x, —t). 

Non- equilibrium dynamics, instead, can be studied by leaving T finite and by 
making the initial distribution Vic explicit [T0| [H]. A typical choice is a Gaussian 
weight in which case can still be absorbed into a redefinition of the fields and g. 
Stationarity is lost out of equilibrium and correlation and linear response functions 
depend on all times involved in their definitions {t and t' in Eqs. (12.101) and (12. lip ). 
Moreover, the FDT is no longer valid [TU [33l 

In addition to R defined in Eq. (I2.10p and C defined in Eq. (12. lip , one can construct 
the quadratic correlator (0j(x, t)(j)j{x', t')) which, independently of the color of the noise, 
vanishes identically due to causality. 

2.1. Scaling 

In the case of stochastic dynamics with white noise, a systematic RG analysis confirms 
the phenomenological scaling behavior of the linear response and correlation functions 
both for T — )■ oo and T finite corresponding, respectively, to equilibrium and non- 
equilibrium relaxational dynamics. In terms of the equilibrium correlation length 
~ |r — rd"*^, where is the critical value of the parameter r in Eq. (12. ip . and 
of a dynamic growing length ^{t) ^ t^/^, one expects [TOj [Ti] 

R{M) = [m/mf fnipum/um/m) (2.13) 

and 

c{p,t,t') = p-'^^ m/mr^'-''^ fdpum/um/m). (2.14) 

for the Fourier transform in space of R{x, t, t') and C{x, t, t'), respectively, with the white 
noise value a = 1 [10]. In the previous expressions, u is the standard static critical 

|j Note that the Ito prescription of the Langevin equation (j2.2|) implies 0(0) = in the stochastic path 
integral description [37l |3] . 
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exponent associated with the correlation length, z is the dynamic critical exponent 
which characterizes the different scaling behavior of space and time, whereas 77 is the 
static anomalous dimension of the field (p cind it controls the power-law spatial decay 
of the static correlation function. 6 is the so-called initial-slip exponent [T0| [T3| [H] 
that accounts for the effects of the initial condition in the case of finite T. It is a 
novel universal quantity if the relaxation occurs from a disordered initial state, whereas 
it is related to known equilibrium exponents if the initial state has a non-vanishing 
average value of the order parameter [TT], [12]. In Eq. fl2.13p and Eq. f l2.14p f^^c are 
scaling functions which become universal after the introduction of proper normalization. 
Equilibrium dynamical scaling is recovered in the limiting case ^{t') ^ ^{t) ^ ,^cq (i-e., 
in the limit of long times t, t' with finite t ~t'), whereas aging phenomena are expected 
to emerge for ^{t),^{t') <^ ^cq and, in particular, right at the critical point r = re- In 
the presence of specific instances of correlated noise we expect the scaling behavior in 
Eq. fl2.13p and Eq. fl2.14p to be modified both as far as the exponents and the scaling 
functions are concerned. The changes appear at the level of the Gaussian theory and 
non-trivial effects survive in the presence of interactions for certain noise correlations, 
as we shall explain in Sees. 3 and 4. 

3. Equilibrium dynamics 

According to the interpretation of the Langevin dynamics in Eq. (12. 2p as resulting from 
the coupling to an equilibrium thermal bath, after a sufficiently long time the system is 
expected to relax towards an equilibrium state characterized by the effective Hamiltonian 
"H, i.e., by the static 0^-theory. This relaxation occurs generically and for arbitrary initial 
conditions as long as the asymptotic values of the control parameters of the system (r in 
the case we are concerned with) imply for neither a spontaneous symmetry breaking 
nor criticality, which would indeed provide instances of aging (see, e.g., [39]). However, 
the existence of a wide region of parameter space (r > 0) for which equilibration occurs, 
allows us to conclude that all static properties of a theory with effective Hamiltonian 'H 
carry over to the dynamic field-theoretical action S [see Eq. (12. Sp and Eq. (12. 9p ] which 
generates the dynamic correlation functions and therefore the static ones as a special 
case. The upper critical dimensionality above which the Gaussian theory becomes 
exact is therefore the same as in the 0^ theory, i.e., Dc = 4: (see, e.g., [37]). Analogously, 
the same applies to the static exponents u and rj. In this section we show how this 
arises within perturbation theory. In particular, we determine the conditions under 
which the critical dynamics is modified by the colored part of the noise with special 
focus on the emergence of a cross-over line ac{D, N) which bounds the region within 
which the dynamic exponent z is affected by the color of the noise. We calculate this 
exponent in the white and colored noise cases. 
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3.1. Gaussian theory 

In the T — )■ oo limit the Gaussian part of the action So can be diagonahzed via a Fourier 
transform of the fields defined in Eq. ( lA.ll) . One obtains 

where we used a vector notation = {(j){p,u) , (f){p,u))'^ for the 2A^-component field 
and we introduced the correlation matrix 

Here and in what follows we denote a function and its Fourier transform with the same 
symbol, the difference being made clear by their arguments. In Eq. (13.21) Fj^ stands 
for the Fourier transform of 0(t)F(t) [the G(t) factor is a consequence of the causal 
structure of Eq. (12. 2p ]. As usual, * denotes the complex conjugate. For the colored 
noise in Eq. (12.61) one finds 

^iuj = 7|w|°~"^ [sin(7ra/2) - i sign(a;) cos(7ra;/2)] + 7^. (3.3) 

[Note that for a > 1 a short-time cut-off has to be introduced in order to transform 
Eq. (12.61) . However, the dynamic properties we are presently interested in are determined 
by the leading behavior at small u, which is not affected by the introduction of such 
a cut-off and is correctly captured by Eq. (13.31) . Accordingly, we shall use this form 
irrespectively of the value of a.] In this expression we have added a supplementary white- 
noise vertex 7^ for reasons that will become clear in the following [note that the cut-off 
that has to be introduced in order to make Eq. (12. 6 p integrable for a > 1 effectively leads 
to this supplementary white- noise vertex]. The propagator^ are deduced by inverting 
C: 

1 

iuTi^ + p2 

and 



Ro{p,uj)Sij = {^i{p,u)^j{-p, -uj)) = „ Sij (3.4) 



Co{p,uj)6ij = {(j)^{p,uj)(f)j{-p,-uj)) 

T + F* 

oo'T^^T*^ + toj{p^ + r)(F,^ - F*J + (p^ + ^ ' ' 

By construction they satisfy the FDT [see Eq. (I2.12p ] that in the frequency domain 
reads: 

2ilmRQ{p,u) = —iuCo{p,u). (3.6) 

We recall that we absorbed the temperature in a redefinition of the fields and the 
coupling constant g, and that Co{p,u) is a real function. 

% In what follows we denote the response and correlation function by R and C, respectively. The 
various propagators and quantities within the Gaussian approximation are denoted by the subscript o- 
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The static correlation function Co{p, t = 0) can be obtained by integrating Eq. fl3.5l) 
over the frequency u and, as expected, the result agrees with the static Gaussian 
correlation that one would infer from the Hamiltonian "H [see, c.f., the calculation leading 
to Eq. (]B.4p ]. Consequently, the static critical exponents u and 77 are not modified at 
this order by the dynamics and they take the Gaussian values z/q = 1/2 and 770 = 0, 
respectively. 

We anticipate here that in Sec. 14. H while discussing the non-equilibrium dynamics 
of the present model, we consider the Laplace transform [see Eq. ( ]A.2|) ] of Eq. (12. 2 p 
with g = and the colored noise given in Eq. (12. 6p (i.e., with 7w = 0). This allows us 
to determine the Laplace transform of the response function i?o, formally obtained by 
replacing iu with A in Eq. (13. 4p : compare Eqs. (lA.ip and f lA.2p . This transform can be 
inverted to a form given in terms of the so-called generalized Mittag-LefHer functions 
Ea,i3 defined in Eq. ( IB.ISP and provides a closed expression for RQ{p,t): 

Ro{p,t) = e(t)^E,,„(-At77), (3.7) 

7 

where A = p"^ + r. The equilibrium correlation function Cq is readily determined from 
this expression via the fluctuation-dissipation theorem Eq. (12.120 (see, c.f., App. IB.2l for 
details) : 

Coip,t) = jE^i-A\tr/^) (3.8) 

where Ea{z) = Ea^i{z). In Fig. [U^a) we plot ACq as a function of the (dimensionless) 
scaling variable u = \t\{A/'yY^°' associated with time t. For a — > 1 one recovers the 
purely exponential dependence e~" (indicated by the decreasing dashed curve in Fig. [1]) 
which characterizes the case of white noise. As a decreases, instead, the correlation 
function displays a faster initial drop followed by a slower decay at large values of 
u. Indeed, taking into account the known asymptotic behavior of the Mittag-LefHer 
functions [c.f., Eq. (1B.19P ]. these curves decay algebraically as ~ l/[r£;(l — a)u°'] for 
u — 00. In panel (b) of Fig. [T]we use a log-log-scale to compare the curves shown in 
panel (a) with their corresponding leading asymptotic algebraic decays, indicated by the 
straight dashed curves for u > 5. As a — the approximation provided by the leading 
term of the asymptotic expansion becomes less accurate in this time span and one needs 
to go to longer times to reach the asymptotic regime. The curves in Fig. [T] clearly 
illustrate the crossover between an exponential and an algebraic asymptotic behavior of 
the correlation function as a decreases below the value a = 1. 

For the generic case of the noise in Eq. (13.31) with 7, 7^ 7^ 0, the propagators Rq and 
Cq do not have a simple analytic form in the time domain, in contrast to the familiar 
exponential relaxation which characterizes the case with white noise (7 = 0, 7„ 7^ 0) 
briefly recalled in Eqs. fIB.ip and flB.2p and to the purely colored case discussed in the 
previous paragraph (7 7^ and 7w = 0). In spite of this difficulty, the Gaussian value 
zq of the dynamic exponent z can be determined by comparing the scaling of the first 
two terms in the denominator of Rq for small u and p since one expects u ^ p^ from 
the definition of z (see, e.g., [37]). First we note that for small u, Eq. (13. 3p scales as 
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Figure 1. Scaling function ACq of the Gaussian correlation Cq in equilibrium 
(T — > oo) as a function of w = |f:|(A/7)^/" for various values of a, with A ^ p^ + r. (a) 
The horizontal dashed line corresponds to the limit a — > 0, whereas the other dashed 
line points out the purely exponential behavior e~", which is recovered for a = 1. 
The solid lines, instead, correspond to a = 0.2, 0.4, 0.6, and 0.8, from bottom to top 
at small u and in the reverse order at large u. (b) Log-log plot of the curves shown 
in panel (a) compared to their corresponding leading asymptotic algebraic behavior 
inferred from (jB.19[) . which are indicated as (straight) dashed lines. 



\u\ 



|"~^ for a < 1, whereas Fj^^ ~ 1 for a > 1: in the former case the effect of 
the colored part of the vertex is dominant, whereas in the latter the contribution of the 
white noise (oc 7w) dominates. As a result, from the scaling uTi^^ ~ one can read the 
Gaussian value Zq of the dynamic exponent: 



.(col) 



2/a for a < 1 



zo=i , , (3.9) 
^ zS") = 2 for a > 1. 

A similar effect is observed in diffusion processes with colored noise, the so-called 
fractional Brownian motion |26]. The particle's displacement is stationary and 
characterized by an a-dependent exponent which is called Hurst exponent in this 
context. 

By rescaling the momentum p and frequency u according top ^ b^^p and u i— t- b~^u 
with b the scaling parameter of the RG flow, one deduces the Gaussian scaling behavior 
of the response and the correlation propagator. We infer from Eqs. (13. 4p and (13. 5p that 

with a similar expression for Co, where the prefactor on the left-hand side (Ihs) is 
replaced by b^^~~^. As anticipated, one can identify two asymptotically scale-invariant 
behaviors (the so-called Gaussian fixed-points in the parameter space) as the Gaussian 



Critical dynamics with correlated noise 



12 



critical point r = is approached. They correspond to P = (7^ = 0,7 7^ 0) for a < 1 
and Pw = (7w 7^ 0,7 = 0) for a > 1, i.e., to the cases in which either the colored or 
the white noise is relevant. The latter reduces to the standard Model A dynamics |37j . 
In order for P and P^ to be fixed points, it is necessary that the corresponding non- 
vanishing coupling strengths, either 7 or 7„, are constant under renormalization which, 
as expected from Eq. (13.91) . implies z = 2;q™^'' = 2/q; for a < 1 (P) and z = z^^ = 2 for 
a > 1 (Pw). 

In order for the action So to be invariant under the momentum and frequency 
rescaling discussed above, one has to rescale the fields (pi and 0j as (f)i{b^^p,b^^"ijj) 
b'^'^(j)i{'p,u) and (j)j^{b~^p,b~''"uj) t-j- b'^'^(j)i{p,uj) where and d-^ are the so-called scaling 

dimensions of the fields and 0, respectively, in the (p, a;)-domain. (Below we shall 
introduce the scaling dimensions of the fields in the time-domain; in order to keep the 
notation as simple as possible we do not include an additional subscript to distinguish 
the two cases but we explain in the text which one we use in each case.) The latter take 
the Gaussian values 

d^,o = iD + 2)/2 + zo, (3.11) 
rf^,o = P + 2)/2. (3.12) 

In the white-noise case zo = 2 we recover the standard scaling dimensions of Model A 
critical dynamics [37]. As far as the transformation properties of the propagators under 
these rescalings are concerned we have b~'^'^'i'^^~^^°CQ{b~^p, b~^°uj\ ■ ■ ■) = Cq{p, u; . . .) and 
5-<^0-"'^+^+^o^^(-5-ip'^ ^-zo,^. . . .) = Rq{p^uj; . . .) where the factor 6^+^° comes from the 
^-function which guarantees the conservation of momenta and frequencies. We have not 
specified the scaling of the parameters r, 7 and 7w to lighten the notation. By comparing 
with the scaling behavior of the Gaussian response in Eq. (I3.10p and of the correlation 
function, one confirms the Gaussian values (13. lip and (I3.12p for the dimensions d^ and 
df, respectively. 

In Eq. (13. 3p we added to the colored-noise vertex associated with Eq. (12. 6p a white- 
noise contribution proportional to 7w for the purpose of highlighting the emergence of 
the two distinct Gaussian fixed points P and P„. As we shall show below such a white- 
noise contribution is anyhow generated under the RG flow as soon as one accounts 
for the effect of non-Gaussian fluctuations (i.e., g ^ 0) on the Gaussian fixed-point 
P = (7 7^ 0, 7w = 0) with colored noise alone. 

3.2. The interaction part 

The interaction part of the action reads 

fdudu'du" d^p dV dV ^ y, , n. 

Oint = I TTTT TFT 7 TFT —:<P[—P — P — P , —(^ — UJ — UJ ) 

J 271 2tt 27r (27r)^ (27r)^ (27r)^ 3!^^ y y ) 

in the frequency and momentum domain. Under the naive scaling with Eqs. (13. 9p . 
(13. lip , and (I3.12p one easily obtains the scaling of the coupling constant: g — >■ b'^~^g. 
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The upper critical dimension is thus Dc = 4: independently of a and the effects of 
fluctuations beyond mean-field can be accounted for by using a standard perturbative 
expansion in terms of e = 4 — 

In the presence of the interaction Sint, the scaling dimension of the fields and 
the coupling constants are altered. In addition, we shall show that the crossover 
value = 1, which separates the colored- noise-dominated case from the white- noise- 
dominated one, acquires a dependance on D, thus dividing the [a, D)-plane (for 
fixed) in two distinct regions, each one with different scaling properties. Under a RG 
flow with scaling parameter b > 1 the noise strengths 7 and 7^ scale as 

7w ^ b^-""'^"^- 7w, (3.14) 

which generalize the corresponding Gaussian scaling behavior of these parameters — 
encoded in Eq. f l3.10p — via the introduction of suitable anomalous dimensions 77^ and 
?7„ of 7 and 7^, respectively. These anomalous dimensions 77^ and ?7„ determine the 
corrections to the Gaussian value zo of the dynamical exponent z and the crossover 
value ttc which separates the different regions in the (a, D)-plane. Indeed, let / be a 
length scale and r be a time scale. Dimensional analysis implies t ~ r and x ~ /. From 
Eq. (13. 4 p we infer that 7 ~ r"//^ and 7„ ~ r//^. Consider the case in which the colored 
noise dominates, which corresponds to having 2 — azo + arj^ > 2 — 20 + '7w in terms of 
the dimensions of the noise strengths [see Eqs. (I3.13P and (I3.14p ] with zq = 2/ a. By 
choosing r" = P'j we have t ~ /^/"7^/". Therefore, under an RG flow with / 1— bl 
{b > 1) we deduce from Eq. ( 13131) that t ~ b^/"+''-< /2/"^i/". Qn the other hand, by 
noting that the dynamic exponent z is defined through t — b^t we can readily identify 
the dynamic exponent 2; = 2/a + 77^ in terms of 77^. In the white-noise-dominated case 
we choose r = 1^'-/^ and a similar argument yields the white-noise result z = 2 + rj^. In 
short, 

'2 + r/„ for a > ac{D,N), 

(3.15) 

2/a-|-?7^ for a<ac{D,N), 

and therefore one needs to calculate rj^ and 77^ in order to determine z. 

In the presence of non-Gaussian fluctuations, the scaling dimensions = d^^^ — 
ZQ — rj/2 and d-^ = d^^ — ZQ — rj/2 in the (p, t)-domain of the fields and 0, respectively, 
differ from their Gaussian values by the corresponding anomalous dimensions rj and r] 
(the extra —zq comes from the conversion of df)^^ and d^-^ from the frequency to the time 
domain). In order to determine the resulting scaling in the (p, a;)-domain one has to 
take into account the integral over time that carries a dimension z (which differs from 
the Gaussian value zq)] therefore 

The FDT implies a relation between rj^, r]^,r] and r/, which allows one to express z in 
terms of the latter two. Indeed, the right-hand side (rhs) and the Ihs of Eq. (I2.12p should 
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have the same scahng dimensions; therefore z = d^ — d-^ in terms of the dimensions of the 
fields in the time-domain. Using now the expressions of the field anomalous dimensions 
provided above, transforming into the dimensions in the frequency domain, and replacing 
the Gaussian values in Eq. (13. lip and Eq. f l3.12p one concludes that 

z = zo^ —■ (3.18) 

3.3. Perturhative expansion 

As we explained above, one does not expect any modification of equal-time correlation 
functions, as they are determined by a static theory with the effective Hamiltonian "H in 
Eq. (12. ip . Hence, we focus on the dynamical exponent 2;, the corrections to which can 
be obtained on the basis of the standard perturbative method consisting in a combined 
expansion in the coupling constant g and in the deviation e = A — D from the upper 
critical dimensionality of the model [381 E3 SOl IH]- In performing such an expansion 
one also takes advantage of the fact that g will eventually be set to its fixed-point 
value g* = 0{e). We remind here that the inverse temperature /3 has been eliminated 
by a suitable redefinition of the fields and the coupling constant g. In the following 
we concentrate on the one-particle irreducible vertex functions E] with n external 
0-lines and n external 0-lines, denotecQ by 

The subscripts indicate the order in the perturbation series. For example, V;^'" includes 
all terms proportional to g^, ge and e^. The Feynman rules of this perturbative expansion 
are those associated with the statistical weight e"*^ in Eq. (12. 7p and they are the 
same as in the white noise case [lOl [37], the only difference being in the form of the 
Gaussian response and correlation functions. In the diagrammatic representation of the 
perturbation series we shall indicate the relevant propagators and vertices as depicted in 
Fig. [2j Note that the noise vertex Ti^^ + T*^ [see Fig. M^d)] is diagonal in frequency space 
(i.e., it amounts to a multiplication by an cj-dependent factor) whereas it is non-local 
in the time domain. In addition, we point out the fact that in principle the correlation 
function can be obtained in the frequency domain as a multiplication of two response 
functions by the noise vertex, which corresponds to a convolution in the time domain. 

3.3.1. Renormalization of the noise vertex. Our interest here is to know whether the 
correlated noise modifies the critical behavior of the model. Within the Gaussian 
approximation z is given by Eq. (13. 9p . where we assumed that a white-noise vertex 
is generated under renormalization, a fact that yields two distinct fixed points P and 
P^: the former is characterized by the colored noise and is stable for a < ttc = 1, whereas 
the latter is characterized by the white noise, is stable for a > ac, and it reduces to the 
standard Model A dynamics. We shall show that, on the one hand, expanding around 

+ Our notation differs from the standard one, that is F"'" for the IPI-vertex functions, in order to 
avoid confusion with the noise kernel. 
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(a) (b) 



9 r,, + 

(c) (d) 



Figure 2. Diagrammatic elements of the perturbation theory: (a) response 
propagator, (b) correlation propagator, (c) interaction vertex and (d) noise vertex. The 
straight parts of each line are associated to fields (j>, whereas wiggled lines correspond 
to (p fields. 



P (with 7w = 0) renormalization indeed generates a supplementary white noise vertex 
7w 7^ and, on the other hand, such a vertex becomes relevant at a D- and iV-dependent 
value ac{D, N), where a^D, N) shows corrections to the Gaussian cross-over occurring 
at ttc = 1 for D < A. 

The first correction to the noise vertex V2'^ is given by the second-order diagram 
depicted in Fig. [3] which can be conveniently written as the Fourier transform of its 



Figure 3. Lowest-order pcrturbativc contribution to the noise vertex. 



expression in the time and space domain 



9\N + 2) 
18 

g\N + 2) 



Vf (g, a) = - '-^^^ I d^x I dt e--^-- dix, t) 



+ 00 



dt cos(o-t) / d^x e'^-"^ Cl{x,t), (3.20) 



9 

where the iV-dependent prefactor accounts for the combinatorics of the graph (see, 
e.g., [3^) and Co is the Gaussian correlation function with 7^ = 0. In the last line 
of this equation we used the symmetry C{x,t) = C{x,—t). Since we are interested 
in the critical dynamics, we set r to its critical value r^ = 0{g) (determined, e.g., by 
the value of r at which C{p = 0,t = 0) diverges [371 SOI E])- However, at the order 
g"^ we are presently interested in, one can neglect the shift of the critical point and 
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set r = 0. The leading behavior of the noise vertex is completely determined by the 
small-g and small-cx asymptotics of V^'^. We can set g = from the outset, while the 
small-frequency limit a — >■ has to be considered with care since the tree-level noise 
vertex Vo'^(cr) = 2Rer,„(7„ = 0) = 27|(t|"-^ sin(7ra/2) [see Eq. (Q] diverges in this 

2 

limit for a < 1. At the end we shall see that no contribution to V2' proportional to 
|ct|"~^ is actually generated. In what follows we only take the limit a — when it 
becomes manifestly possible. In this formulation, divergencies arise due to the singular 
behavior of Co at small distances and times, |x| — ?■ and t — t- 0. In order to regulize 
the theory, we introduce a short-distance cut-off i, below which the description in the 
continuum is no longer considered to be realistic. For example, in lattice models, the 
cut-off £ is naturally identified with the lattice spacing. Analogously, a short-time cut- 
off is introduced in the convenient form i^, which is motivated by the scaling behavior 
discussed above. By using, c.f., the scaling form ( IB. 51) of the Gaussian correlation 
Co influenced by the colored noise (see App. B), the second-order contribution to the 
regularized vertex function takes the form 

(3.21) 

where Ad = 271^/"^ /V e{D /2) is the solid angle in D dimensions. 

The Wilsonian renormalization scheme (see, e.g., |10]) amounts to a resummation 
of the perturbation series which is performed according to the following steps: 

(I) Effective vertex functions for the 'slow' fluctuations are determined by 
performing an integration (averaging) over 'fast' fluctuations, within a spatial shell 
between £ and M and at a temporal scale between £^ and {MY . As a result of this 
integration the effective vertex functions — and therefore the coupling constants which 
characterize them — acquire a dependance on the scaling parameter fe > 1. To be more 
specific, consider the typical integral which arises in loop calculations and which can be 
written in the generic form 



/•oo /"OO 

X(£) = / dt I d^X J^(f,t), 



with some integrand J^{x,t). The contribution of the integration over the fast 
fluctuations is then equivalent to X(£) — I{M), an expression which we shall use 
repeatedly below. In the limit 6 — 1 with 6 > 1 one has X{i) — X{b€) = 
- [dl{i) /d In i]\nb + 0{\n^ b) . 

(II) The effective vertex functions calculated in step (I) depend on a new cut-off M. 
In order to recover the original cut-off i one rescales the coordinates and fields in the 
frequency and momentum domain according to 

T ^ 5D/2+l+.~.o-rJ/2T 
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The resulting expression is multiplied by h~^~^ which accounts for the rescaling of the 
integration measure in the Hamiltonian. 

(Ill) In order to study the evolution of the coupling constants under the 
renormalization procedure it is convenient to consider the case h ^ 1^ which corresponds 
to an infinitely thin integration shell. In this case the evolution equations for the coupling 
constants are coupled differential equations that depend upon a and the anomalous 
dimensions introduced by the rescaling in step (II). The anomalous dimensions are 
determined by requiring that all coupling constants have a finite asymptotic value under 
the RG transformation for 6 — oo. 

Applying step (I) to the noise vertex function V'^''^ we derive Eq. f l3.2ip with respect 
to In I and we multiply the result by In h. By defining [see Eq. fIC.ip 

(3.3, 

with 

u = Az)(?/(27r)^, (3.24) 

we find that the effective noise vertex V'^'^(0, a — )■ 0; 6^) for the slow fluctuations with 
short-time and -distance cut-offs hi and (fe^)^, respectively, is given by 

V°''(0, a^Q-W)= - (r..^o + a^o) - «'^°''(0; 7) In 6 

+ 0{u^\Yi^h,u^). (3.25) 

For details on the calculation of £^°'^((t; 7) we refer to App. O 

Clearly, the form of the effective noise vertex has changed, as the term £^°'^(0;7) 
generated by the non-Gaussian fluctuations has the form of a white-noise contribution, 
whereas the coefficient 7 of the colored noise is not modified up to this order in 
perturbation theory. As a result, it is convenient to account for the contribution of 
a white- noise vertex from the outset, by replacing Fjo- by Fj^- + 7^. This implies that 
the Gaussian correlation function Cq that determines the loop correction still has a 
scaling form but with a scaling function gc^ that is now a function of two variables, see 
Eq. f lB.12p . The correction that is generated depends on both 7 and 7„, we denote 
it by £^°'^(0; 7, 7w) and we explicitly calculate it in Eq. (IC.ll) . 

The effective noise vertex depends on the cut-off M. Following step (II) of the 
renormalization procedure we rescale the effective noise vertex as specified in Eq. ( I3.22p . 
The coupling strengths of the colored and the white noise 7 and 7^ become running 
coupling constants 7(6) and 7w(&) and in the limit 6 — 1 they satisfy the set of coupled 
differential equations 



^7 
(9 In 6 

and 

^7w 
(9 In 6 



azo - -{rj - 7]) - 7] 7 + C(e') (3.26) 



2 

f] + r] 



7w + V2£:O'^(O;7,7w) + 0(e^), (3.27) 
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valid at the critical point, u* = 0{e) is the fixed point value of the coupling constant, 
i.e., the value at which the effective coupling constant u{b) — obtained by applying 
the procedure outlined here to the 4-point function — flows for 6 — )■ oo and D < 4. 
For D > 4, u* = and the scenario within the Gaussian approximation presented in 
Sec. 13. H is not altered by the interaction. Accordingly we focus below on the case D < 4. 
Two additional differential equations can be written by considering how the coupling 
constant u in V^''^ and the coefficient of the term cx: in V^'^(g, . . .) are modified by the 
non-Gaussian fiuctuations. In particular, the requirement of an effective 6-independent 
coefficient of fixes rj to its well-known static value [37] (see Sec. 13.3.31 for further 
details) . 

In order to determine the critical exponents we demand that the amplitude of the 
noise vertex in the effective Hamiltonian be constant as explained in step (III) of the 
renormalization procedure. Neglecting for a while the contribution of the non-Gaussian 
fluctuations to Eq. fl3.26p and Eq. f l3.27p (which amounts to setting u* and the anomalous 
dimensions to zero), one can easily solve them and recover the Gaussian picture which 
we anticipated in Sec. 13.11 Indeed, 7(6) ~ whereas 7„(&) ~ fe^-^o g^g ^ which 

implies 'y(b)/'j^(b) ~ b^'^~°'^^°. Independently of the value of zq > 0, this ratio tends to 
zero for a > 1. The associated flxed point is characterized by a flnite 7w(^ — ^ 00) with 
a vanishing 7(6 — )• 00), i.e., the flxed point introduced in Sec. 13.11 In order for 7w(&) 
to stay flnite, it is necessary to have zq = z^^ = 2 in Eq. f l3.27p . as expected from our 
previous discussion. On the contrary, for a < 1, 7(6)/7w(6) — )■ 00 for 6 — )■ cxo and the 
associated flxed point has a flnite 7(6 — 00) and a vanishing 7„(6 — > 00), corresponding 
to the flxed point P of Sec. 13.11 The former condition requires zq = ^g™''* — 2/a in 
Eq. (I3.26p . consistently with the discussion therein. 

Including now the effects of non-Gaussian fluctuations, the colored-noise flxed point 
P with Zq = zl^"^^ = 2/a and 7(6 — ?■ 00) 7^ is characterized by a value of rJ such that 
the Ihs of Eq. f l3.26p vanishes. This yields 

, . . (1 - ^) „ . ^ (1 - f ) 5|±|^.^ . . (3.28, 

We replaced 77 by its static value given in [371 SD] since, as we shall show in Sec. I3.3.3[ 
it is a-independent. The value z^^°^^ of z at this flxed point is determined via Eq. ( ]3.18p 

' a a[ 4(iV + 8)2 

The flxed point P is stable in the (a, D)-plane (region C in Fig. H]) as long as the value of 
7w(fc) determined by Eq. f l3.27p at the flxed-point P with rj = and zq = Zq™^-* = 2/a 
stays flnite for b — )■ 00. A crossover towards the flxed-point Pw in the (a, D)-plane 
(region W in Fig. H]) controlled by the white noise occurs as soon as j^{b — t- 00) — 00. 
In this limit, £^°'^(0;7,7w — > 00) ~ 7w^^^'^ independently of 7 [as long as 7(6) remains 
finite, see Eq. ( IC.SP for details] where the constant S^''^ is given in Eq. (1C.6P and is such 
that 

«*^^r = ^^^31nie^ + 0(a. (3.30) 



+ Oie^). (3.29) 
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Thus, the equation which determines the evolution of 7„ at the fixed point P becomes 



^7w 
d\nb 



2 - - 



^(col) 



r] z 
— + 



(col) 



*2(70,2 



-U 6, 



7^^ 



(3.31) 



and the crossover occurs as soon as the the quantity in brackets changes sign. The 
expression of the crossover fine is readily determined by taking into account the values 
of 4™'\ and reported in Eqs. ^JM) . ^JM . and (l330|) : 

e' + C(e=^). (3.32) 



2 3 



(A^ + 8)2 

For a > ac (region W in Fig. Hj), 7„(6 — )■ oo) — )■ cxo and the point P is no longer 
a fixed point of the rescaled effective action, as the white-noise contribution becomes 
predominant. In order for it to become constant and therefore to determine the fixed 
point Pw, ^ in Eq. f l3.27p should now take the value rf^^ such that d^^jdXv^h = 0, 
with zq = 2q"^ = 2. Assuming that the coefficient 7(6) of the colored noise vanishes 
asymptotically for 6 ^ 00, the Ihs of Eq. (Km becomes -(r/(™) +r])/2 + {z/2)u*^S^''^, 



where we used the fact that £^°'^(0;7 
the same equation vanishes implies 



0,7. 



7wi^w^- The condition that the rhs of 



Tj = Tj^ ' = —rj 



2n*^£ 



0,2 



N + 2 
(iV + 8)' 



6 In 

3 2 



and from Eq. fl3.18p . 



2 + 



N + 2 



, 4 1 

3 n 

3 2 



(3.33) 



(3.34) 



in agreement with [3l[37]. In order to verify the consistency of the assumption 7(6) — 
for b ^ 00 under which Eq. fl3.33p has been derived, one can specialize Eq. f l3.26p to 



the white- noise fixed point P^, by using rj 



^(w) (w) 



2, and z^^^ [see Eq. f l334D ] as the 



values of f/, zq, and z. Accordingly, the term in parenthesis in the rhs can be written as 
—2{a — ac) + 0{e^) and therefore 7(6) ~ b'"^^"'""^ indeed vanishes for a > ac as 6 — )■ 00. 
This also proves that the white-noise fixed point Pw is stable against the perturbation of 
the colored noise as long as a > ac, a statement which complements the one presented 
above about the stability of P. 

Summarizing, Eq. fl3.32p determines the line in the {a, il')-plane which separates 
region W from region C: in the former, the white noise dominates and z = (in 
agreement with [21 EZ]), whereas in the latter the colored noise dominates and z = 
is given by Eq. fl3.29p . Figure H] illustrates this scenario for = 1, 4, 00. 



3.3.2. Renormalization of the self- energy and FDT. The fiuctuation-dissipation 
theorem (FDT) expressed in Eq. 02.121) is a consequence of a symmetry of the action in 
equilibrium [TOl [35] and it has to be preserved under renormalization. Therefore, the 
noise vertex and the memory kernel have to be related by the FDT even beyond the 
Gaussian approximation, that we analyzed in Sec. 13.11 Here we explicitly show that 
this relation is still valid when non-Gaussian corrections up to 0{e^) (or, equivalently. 
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Figure 4. Boundary between the regions W and C of the {a, Z))-plane characterized, 
respectively, by white and colored noise. The boundary curve a = ac{D,N) as a 
function of the spatial dimensionality D is reported here for = 1 (solid line, Ising 
universality class), 4 (dashed), and oo (dotted) where the ©((4 — Z?)'^)-correction in 
the corresponding perturbative expression p.32p for D < 4 has been neglected. The 
vertical dashed line indicates the lower critical dimensionality of the model for iV > 1 . 
The coefficient of the term 0{{4—D)^) in Eq. p.32p takes its maximum value for = 4 
(dashed curve) and then it decreases monotonically as a function of A^, vanishing for 
N ^ oo. For D > A, ac takes the D-independent Gaussian value a^o = 1 (dotted 
line). Clearly, the dependence of the boundary curve on the dimensionality D is 
quantitatively rather weak. 




Figure 5. Second-order contribution to the self-energy. 

0{g'^)) are accounted for. The first correction to the memory kernel comes from the 
second-order self-energy contribution 

Vl'\q, a) = j^^t j d^x e"^-'^-'''' Cl{x, t)Ro{x, t) (3.35) 

represented in Fig. Note that Rq is causal and restricts the time integral to run 
over positive values only. The expansion of this expression as a power series in a and q 
allows one to identify the terms which contribute to the renormalization of the different 
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parameters of the Gaussian vertex VQ'^{q,a) = icrTifj + + r. The terms which are 
independent of both a and q contribute to the renormahzation of the parameter r (which 
is also modified by an 0{g) term not discussed here), the terms proportional to a^q"^ 
contribute to the renormahzation of the fields and those proportional to ia^q^ to the 
renormahzation of the memory kernel F. First of all we observe that the FDT given 
in Eq. fl2.12p allows us to express Rq in Eq. f l3.35p as —dfCo. An integration by parts 
yields 

V'2\<1,^) = -^^^^ J e^'^-'!^C',{x,0)-za J^^dte-^'^'CUx,t)Y 

(3.36) 



Hence, 



18 J Jo 



a 



^vriq,^), (3.37) 

where the last equality follows from a comparison with Eq. f l3.20p and shows that the 
FDT in the frequency domain [see Eq. f l3.6p ] is satisfied by the corrections 0{g^). 
Note that the vertex V^'^ receives also a correction of 0{g) given by a tadpole 
diagram which, however, is a real constant and does not contribute to the imaginary 
part. We conclude that up to and including the second order in the coupling constant 
2 Im V^'^(0, cr) = —(xV^'^{0, a). (This proof can be readily extended to the corresponding 
regularized vertex functions, characterized by short- time and distance cut-offs.) 



3.3.3. Renormalization of the self-energy: the anomalous exponent rj. In the same 
spirit as before we can deduce the first correction to the static exponent 77. It is 
instructive to see why the dependence upon a does not affect the final result, even 
though the regularized expression of Eq. fl3.35p does via Cq and Rq. In order to single 
out the contribution of Vg'^ to the coefficient of g^, one expands Eq. fl3.35p — suitably 
regularized as discussed above — up to second order in g, finding 

Vl'\q,a = 0-J) = -q'g'^—^ dt dx ^^(x , t)/2o(^, t) + • • • , (3.38) 

where the dots indicate all the terms which do not contribute to the field renormalization, 
i.e., which are not proportional to a^q'^. In Eq. f l3.38p we used the fact that, for a 
generic function /, J d^x XiXj f{\x\) = {6ij/D) J d^ x {xl"^ f {\x\) , which is valid also for 
the regularized integral. As in the case of Eq. fl3.36p one can take advantage of the 
FDT, Eq. f l2.12p . to express the integrand in Eq. f l3.38p as a total derivative, which can 
be integrated to yield 

V'Aq. 0; i) = qV ^^4^ r "^"^ ^° + • ■ ■ ^^-^^^ 
We note here that even though the (full) dynamic correlation function C (and therefore 
its Gaussian approximation Cq) depends on the value of a, the static correlation function 



Critical dynamics with correlated noise 



22 



C{x,t = 0) does not. This is explicitly shown for Co in Eq. (lB.4p . While the limit 
£ — 7- of the rhs of Eq. fl3.39p cannot be explicitly taken due to the short-distance 
singularity of the integrand, such a limit can be taken for the correlation function, i.e., 
Co(x, t = i^) ~ Co(x, t = 0) and therefore the expression of V2^{q,0;i — )■ 0) becomes 
— as expected — independent of a at the leading relevant order in i. Applying the 
same renormalization procedure as in Sec. 13.3. ![ Eq. f l3.39p can be used to calculate the 
effective vertex V2^{q, 0; bi) after having integrated out the fast fluctuations. Similarly 
to Eq. fl3.23p one defines 

,V." + ,„.-M;fM, (3.40) 

For 6 — )■ 1, the resulting effective vertex is 

Vl'\q, 0; M) = q^ + q\^S^'^ In b + 0{u^ In^ b,u^) + ... (3.41) 

In order to recover the original cut-off i we rescale the fields and coupling constants 
according to Eq. fl3.22p and we take 6 — )■ 1. The part of V2^{q, 0; i) that is proportional 
to (T°g^ satisfies the evolution equation 



dlnb 



-q^ [r] - u*^S''^] + 0{u^ In 6, u^) + . . . . (3.42) 



By demanding that the amplitude of V"^^ be constant and by using the numerical value 
of E^'^ calculated in Eq. ( IC.SP we find 

^ = n-^^'^ = ^(^^e^ + C^(a, (3.43) 

i.e., ?7 has the same a-independent value as in the static theory confirming our 
expectations alluded to at the beginning of Sec. [31 



4. Non-equilibrium dynamics 

J^.l. Preliminary remarks 

In this Section we investigate the non-equilibrium dynamics assuming that the model 
is prepared in some initial condition at time t = and that it is let relax subsequently 
at its critical point. This problem has been studied in detail in the white-noise 
case [TOl [T3] . The analysis reveals the emergence of an interesting scaling behavior of 
two-time quantities, usually referred to as aging (see, in this context, [HIISSIES])- More 
precisely, the relaxation is studied via the field-theoretical action S in Eqs. (12. 8p and 
(12. 9p with T = 0, supplemented by a suitable distribution Vic for the initial condition 
at time t = T = 0. In particular, a high-temperature disordered state is modeled by a 
Gaussian distribution with zero mean: 

InVic = - /d^'a; y0'(x,t = 0) (4.1) 

The parameter tq sets the inverse width of the initial distribution. Within the Gaussian 
approximation, the field has a scaling dimension d^^o given by Eq. (13. 9 p in momentum 
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and frequency space, i.e., a dimension d^^ — zq — D in the space and time domain. Using 
this dimension for the initial field (f){x, t = 0) we find that tq — )■ B'^tq under rescaling. 
Consequently, the width of the initial distribution shrinks to zero as 6 — )■ oo, leading 
to a zero effective value of the initial order parameter t = 0) = and, therefore, 
to a correlation function with Dirichlet boundary conditions at t = 0. In App. |D] 
[see Eq. f ID.Sp and Eq. flB.131) ] we show that the Laplace transforms of the Gaussian 
propagators are: 

(A + K,){XTx + + r) 

r + r 

Co{p; A, k) = — — \ ^ (4.3) 

(A + K) [XTx + + r) {kT^, + + r) 

where the Laplace transformed noise is Fa = 7A""^ + 7w with A G 71+ . [In order to 
transform Eq. (12.60 for a > 1 one should introduce a short-time cut-off. However, as we 
pointed out after Eq. f l3.3p . this modification is not necessary as long as one is interested 
in the leading long-time, near critical dynamic behavior of the system. Accordingly, we 
shall use this form for irrespective of the value of a.] As in the equilibrium case, the 
propagators have a simple analytic form in the time domain only for a = 1 or 7„ = 0. 
It is easy to show that the response propagator is the same in and out of equilibrium 
and, therefore, that it is time-translationally invariant. 

The correlation function Cq can always be written as the sum of the Gaussian 
equilibrium correlation C^^^ and the remaining non-equilibrium contribution, which we 
denote by Cq°°^ and which will play an important role in fixing the genuinely non- 
equilibrium properties of the relaxation, e.g., the non-equilibrium exponent 6 and the 
effective temperature. The Laplace transform Cq^^ (p. A) of the equilibrium correlation 
function c'^\'p,t) can be obtained from Eq. ( IB. 211) : 

C^o\v. A) = ^\ l\ . (4.4) 

+ r XTx + p^ + r 

The full non-equilibrium correlator Eq. ( 14. 3 p can be expressed as 

Coip, A, K) = Cl^'\p^>^) + Ct\p.-) _ + ^^ct\p, X)Ct\p, K) (4.5) 

A + K 

which displays the fact that Cq is the sum of an equilibrium time-translationally 
invariant term and a non-stationary term. Indeed, the Laplace transform C[F]{X, n) 
with respect to both t and t' of any (translationally invariant) function F{\t — t'\) is 
given by C[F]{X, n) = {Fx + Ft^)/{X + k), which is exactly the form of the first term 
in Eq. ( 14. 5 p [28] . Accordingly, we can identify the non-equilibrium part Cq'^'^^ of Co as 
Cq^°\p, X, k) = —{p^ + r)CQ\p, X)Cq\p, k) which translates, by virtue of Eq. (]B.2ip . 
into the non- stationary expression [t, t' > 0) 

Ct\p;t,t') = —l^E^{-{p' + r)t'^h)E^{-{p' + r)t'-h), (4.6) 

_|_ ^ 

where we wrote the equilibrium Gaussian correlation function in terms of the Mittag- 
Leffler function, anticipated in Eq. (13. 8p and discussed in App. IB. 21 
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4-2. General non- equilibrium renormalization group analysis 

The addition of the initial condition to the action modifies the scahng of the fields at 
the boundary t = compared to the one in the 'time bulk' t > [lO]. In addition, to 
the bulk renormalization a new initial time renormalization is required, which gives rise 
to contributions 'located' at the time surface (that is the hyperplane determined by the 
condition t = 0). These can be absorbed by introducing a new anomalous dimension of 
the initial field 0o(p) = 'Pi'Py 0) (see |12] for an application to surface critical phenomena). 
The general scaling of the initial fields in the time and momentum domain reads [cf. 
Eqs. f l336|) and (KIT} ] 

^(b-'p^t = 0) ^ 6^/2+i-''/2 0(p,O), _ 

0(6~V,^ = O) ^ 6^/2+i--o-W2-W2 0(^-^0), ^ 
where is a new exponent, with a Gaussian value r/jno = 0. Note that the anomalous 
dimension of the initial response field 0(p, 0) is allowed to differ by r/j^ from its bulk 
value. In [101 E] oiig can find a careful analysis for the white-noise case where it is 
explained why only the initial response field has to be renormalized. Here, we make the 
same assumption and we check its validity a posteriori, r/j^ is related to the so-called 
initial-slip exponent 9 [introduced at the end of Sec. [21 see Eqs. fl2.13p and fl2.14p ] by [10] 

9 = -r]Ji2z). (4.8) 

The analysis in [TOl [I4] has to be slightly modified to deal with colored noise. Our 
starting point is the general leading scaling behavior of the critical correlation functions 
gn,n,no bulk ficlds (f), n bulk response fields (p ^^id uq initial response fields </)g, 
evaluated at the set of points {p, t} in momentum and time: 

where 5{n,n,no) = -D + n{D/2 + 1 - r]/2) + n{D/2 + 1 - zq -r]/2) + no{D/2 + 1 - zq- 
r//2— r/ijj/2). [In writing Eq. f l4.9p and the analogous relations presented below, we always 
understand that the correlation functions on the Ihs and rhs are characterized by the 
different length cut-offs hi and i, respectively.] This scaling behavior is a consequence 
of the scaling dimensions of the fields (f), (j) and (pQ as functions of time and momentum 
[compare to Eqs. fl3.16p . f l3.17p and (14. 7p ] and the dimension of the 5-function ensuring 
the total momentum conservation. Note that all correlation functions with an initial 
field 00 = 0(P) 0) vanish [see the discussion at the end of Sec. 14. Ij . Specifically, the two 
point correlation and response functions (^^'"^'^ and respectively) scale as 

C(p; t, t') ~ C{p/h- bH, bH'), (4.10) 

By choosing b = (t — t')~^^^ these scaling forms become 

C{p-t,t') ^ {t-tT-"^)/^ Fc{{t-t'f'^P.t'/t), 

R{p; t, t') ~ (t - t')(2— Fn{{t- t'y/'p, t'/t). ^ ■ ' 

In general the scaling functions Fc and are not expected to have a finite, non- 
vanishing value for t' — )■ 0. In order to deduce their behavior for small t' we employ a 
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short-distance expansion [57] of the fields 0(p, t') and (j){p, t') around t' = 0. However, 
these are not independent. Indeed, the full correlation and linear response functions 
verify the equations [10] 

C{p;t,t') = J^dsJ^ds' Roip;t,s)V^'\p,s,s')Co{p;t',s') 

+ /o*' ds' f^' ds Co(p; t, s)V^'\p, s', s)i?o(p; t', s') (4.13) 
+ ds /J' ds' Roip; t, s)V0'2(p, s, s')i?o(p; t', s') 

and 

R{p;t,t')= [ ds [ ds' Ro{p;t,s)V^'\p,s,s')Ro{p;s',t'), (4.14) 

where V^'" are the (not necessarily one particle irreducible in contrast to the ones 
introduced in Sec. 13. 3p vertex functions with n amputated external field and n 
amputated external response field legs, respectively. In writing these expressions we 
accounted for the causality of V^'^{p, s, s') oc 6(s — s'). After taking the derivative 
of Eq. f l4.13p with respect to t' only the first term in the rhs survives in the limit 
t' — )■ [note that Rq{p]s,s) = [35]]. By comparing the resulting expression with 
the rhs of Eq. ( 14.14P one notices that the equations differ only by the last factor 
in their integrands, dt'Co{p;t' , s') and Ro{p; s' ,t'), respectively. We deduce that if a 
relation between the dimensions of the time-derivative of the initial field and the initial 
response field exists within the Gaussian approximation, it should be preserved when 
non-Gaussian fluctuations are accounted for. Let us then examine the propagators. We 
focus on region C where they satisfy the equation, 

dt'Co{p;t,t' ^0) ^t'"^^ [ dsT{t- s)Ro{p;s,t' ^0) (4.15) 

Jo 

[proven in App. [Dl see Eq. (ID. 81) ]. In the early t' limit we formally expand the fields 
according to 

(j){p,t' ^ 0) ^ v{t')MP) and ^{p,t' ^0)^Tp{t')Mp}. (4.16) 

is proportional to (f)o{p) and is proportional to 0o(p) since the former vanishes while 
the latter is allowed to be finite for t' — t- 0. We see from Eq. (14. 150 that, under the 
rescaling t — b^t, s — )■ b^s and p — )■ p/b (leaving t' unchanged), the scaling dimensions 
d^^ and d-^^ of 0o and 0o, respectively, verify 

^^0=^(1-") + ^00- (4.17) 
For a = 1 this reduces to the relation found in [10]. The expansion of in Eq. fl4.16p 
can be used to calculate the correlation function Q^'^''^{p;t,t' — )■ 0) ~ '(p(t')Q^''^'^{p;t) 
and by matching the scaling dimensions of the Ihs and rhs with the help of Eq. (14. 9p 
we conclude that (p{t') ~ t'~^ where 6 is given by Eq. (14. 8p . Besides, the rescaling of t' 
(keeping t and s unchanged) implies f{t') ~ if the scaling dimensions of the Ihs 
and rhs in Eq. (14.150 are to match. Hence, the small-t' limit of the response function is 

R{p-t,t' ^0) ^W){^{-p,t)MP)) r^t'-'g''''\p,t), (4.18) 
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where we introduced a short-hand notation for the arguments of Q^'^'^ in which we 
only write the non- vanishing time t. The scahng properties of Q^'^'^{p,t) are given by 
Eq. (HH): 

g'''''\p,t) ^ t-'^v/2+v/2+rj,j2+zo-2)/z gl,0,l ^^1/ z 

= t^2-v-z)/z+e glAl^^l/y^ 1) (419) 

where we used the relation between anomalous and dynamic exponents, Eq. (13.181) . and 
the relation between 9 and r]^^, Eq. (14. 8p . Consequently, taking Eqs. (I4.12p . (I4.17P and 
(I4.19P into account, we conclude that 

R{p- 1, t' 0) ^ Fnit'/y, 0). (4.20) 

A similar analysis of the scaling behavior of the correlation, taking into account 
Eq. dHZD, yields 

C{p; t, t' ^ 0) ~ Fc{t'/y, 0). (4.21) 

These results are used to capture the singular behavior of the scaling functions in 
Eq. (I4.12P by writing: 



(4.22) 



R{p-t,t') ^ {t~t'Y'-^-^'y^ (^L^\^^^t-t'Y/y,t'/t), 
C{p;t,t') ^ {t-t'f-^^/^ (^X \c{{t-ff^p,t'/t), (4.23) 



with 

f 1 for a > aJD,N), 

la for a<ac{D,N), 

which encompass the white noise result a = 1 [TU] for a > ac- The scaling functions 
Fc and are regular for t' — > and depend on a. Moreover, in the RG sense they 
are universal functions up to an overall amplitude and the normalization of their first 
argument. 

The emergence of a 7^ 1 for colored noise can be checked within the Gaussian 
approximation by looking at the initial-slip behavior of the propagators Rq and Cq 
with a < 1. First of all, note that 6 takes the value ^0 = within the Gaussian 
theory, as one can infer by comparing the scaling form Eq. (I4.22p with the expression 
for the non-equilibrium response Rq at criticality, which coincides with the equilibrium 
one in Eq. (13. 7p and is invariant under time translations. Using this value 6*0 of 6 
one has lim kCo{p; \, k) ~ and lim kRo{p; \, k) ~ from Eqs. (14. 2 p and (14. 3p . 
respectively. 



4.2.1. The initial-slip exponent 6. Out of equilibrium the first correction to the 
self energy leads to a modification of the scaling of the initial response field. The 
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response function up to first order in the perturbative expansion reads, for zero external 
momentum, 

R{0;t,t';i) = Ro{6;t,t') + [ ds RoiO;t, s)Be-iis)Ro{0; s,t'). (4.25) 

Jt' 

Bi-i (s) stands for the 'tadpole' diagram represented in Fig. [6], which can be calculated 
by using standard Feynman rules in the time domain [H El [7j, whereas is the large- 
momentum cut-off introduced in order to regularize the otherwise divergent integral 
defining Bi-i{s): 

B,Ms) = -?^f Chip's, s). (4.26) 

6 J\p\<e-^ (2vr)^ 




Figure 6. First order contribution to the non-equilibrium self-energy. 



(ne) 



The renormalization of the initial response field is due to the non-equilibrium part Cq 
of Cq. Indeed, the equilibrium part Cq^^ is characterized by time-translation invariance 
and therefore it contributes with a time-independent function of p to Cq{p]s,s) in 
Eq. f l4.26p . In turn, such a function results in a time- independent contribution B^'^Ji 
to Bi-i, which can be thought of as due to a shift r i— )■ r — of the mass r in the 
expression of the response functions Rq appearing in the rhs of Eq. f l4.25p . i.e., as a 
mass renormalization. [We recall that Ro{p] t, t') actually depends on the two times via 
t — t'.] One can check that this term yields the correct first order correction to the 
critical exponent u which is the same as in the static theory. 

In view of the renormalization procedure outlined in Sec. 13.3. ll we need to calculate 

r%-.BeMt) = - ^(^ + ^) r^ ct\\p\ = r\t,t) (4.27) 

o 

in the limit £ — > 0, for r = and D = 4. By using the asymptotic expansion of the 
generalized Mittag-Leffler functions Eq. fIB.lQp and their definition, Eq. (IB.lSp . one finds 

\ (—xy^ /r e(1 — <y) for X — )• — cx) 

' \^a(0) = 1 for |x| < 1 ^ ' 

and therefore, using Eq. (14. 6p . 

t Go (£ ^ ^ '^^^^^-\or') for rtV7«l. ^ ^ ^ 

Accordingly, £"^c?£-ii?^-i(t) — )■ in the limit £ — )■ for every fixed t > 0. The 
physical interpretation of this fact is that only the initial field is renormalized by 
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Eq. (14.25 p . Indeed the rhs of Eq. fl4.29p for finite i provides an approximation of the 
delta distribution restricted to z E M^, usually denoted by S^{z): 

d{a) 



-6+{z) for i^O, 



where the normalization constant d{a) is given by 



d{a) 



dzEii-z) 



(4.30) 



(4.31) 



and the additional factor 1/2 on the rhs of Eq. fl4.30p has been introduced for later 
convenience in order to have d{l) = 1. Taking advantage of the closed- form expressions 
of the Mittag-Leffler function for a = 1, 1/2, and 0, i.e., Ei{—z) = exp{—z), 
Ei/2{—z) = {2/y/Ti) f^dte^^~^^ , and Eo{—z) = 1/(1 + z) [13], respectively, it is possible 
to calculate the corresponding values of the a-dependent constant d{a). One finds 
d{l) = 1, d{0) = 2 and, after some algebra, d{l/2) = ^/2/n\n{3 + 2^2) = 1.406. . .. 
Hence, 

dBi-i(t) u-i{N + 2)d{a) 



(5+(t°) for £^0. 



(4.32) 



(91n£-i 12 

Using this expression and the one of the zero- momentum response function -Ro(0; t, s) = 
{t — s)"~^/[7r£;(a)] at criticality which follows from Eq. f lB.17p . the derivative of the 
tadpole contribution to the rhs of Eq. fl4.25p can be written as 



d 



ds i?o(0;t, s)Bi^i{s)Roi6;s,t') 



u{N + 2) d{a) 



■Ro{0;t,t') 



(4.33) 



where Sffl = 1 for t' = and otherwise, illustrating the fact that only the 
initial field is renormalized. In deriving this last equation we used the fact that 
5+(r) = 5+(t)/(ar^^). Altogether, the effective response function with cut-off bi reads 



R{6;t,t';b£) = RoiO;t,t') 



uiN + 2)d{a] ^ 
1 + dt'o — , _ _ , , — In b 



(4.34) 



12arE{<y) 

In order to recover the original cut-off i we make use of the scaling relation Eq. f l4.9p 
with t' = 0. By taking into account that t] = r] = 0{e^) we have 

R{6; t, 0; bi) ^ r2+"o+^in/2i?(0; bH, 0; £) 

u{N + 2)d{a) 



i?o(0;t,0) b"-^-/^ 



1 + ~'" ' In 6 

12aTE{o;) 

u{N + 2)d{a] 



(4.35) 
(4.36) 



= RoiO;t,0) 1 + f ln6+ "\-;77 ln6 . (4.37) 
2 12ai E[ct) 

By requiring that the amplitude of the response function be constant at the fixed point 
u* we obtain 



{N + 2)d{a) 



-e + Oie^) 



{N + 8)aTE{a) 

whence we find the a-dependent initial slip exponent from Eq. 



e 



a_ _ {N + 2)d{a) 
'4^'° ~ 4(iVT8jr^ 



(4.38) 



(4.39) 
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Figure 7. Ratio 9/da=i between the initial-slip exponent 9 in Eq. (|4.39p and its 
white-noise value 9 function of a, within the relevant range < a < 1 

at the first order in the e-expansion. The dashed horizontal line indicates the value 
d(l/2) / = ( V2/7r) ln(3+2\/2) = 0.793 .. . corresponding to a = 1/2 (vertical dashed 
line). The dash-dotted line points out the linear behavior 9/9a=i ~ 2a expected for 



In the white noise case a = 1 we obtain 6 = {N + 2)e/[A{N + 8)] -|- O(e^) in agreement 
with the first order result reported in [10]. The dependence of ^ on a is shown in Fig. [71 
6 increases monotonically from ^ = Oata = Oto6' = 6'(l)ata = l which is the 
cross-over value up to 0{e). 

4-2.2. Fluctuation- dissipation ratio and effective temperature. A system which 
equilibrates after a certain finite relaxation time satisfies the FDT. More generally, 
one defines the fluctuation- dissipation ratio (FDR) by 



where is the temperature of the thermal bath (set to 1 in the previous analysis). 
In glassy and weakly driven macroscopic systems with slow dynamics — small 
entropy production limit — this ratio approaches a constant on asymptotic two-time 
regimes in which, moreover, it is independent of the observable used to define the 
correlation and associated linear response and admits the interpretation of an effective 
temperature [301 EI]- For systems with critical points, the asymptotic value 



has been suggested to behave as a universal property [32] and, moreover, as an effective 

temperature, 




(4.40) 




t'— s>oo t— s>oo 



(4.41) 




lOO 



(5X' 



oo 



(4.42) 



Critical dynamics with correlated noise 



30 




Figure 8. Fluctuation-dissipation ratio for the global order parameter (corresponding 
to p = 0) at criticality r = within the Gaussian approximation, as a function of 
the ratio < t' /t < 1 for various values of a. The straight horizontal and diagonal 
dashed lines correspond to a = 1 and a = 0, respectively. The solid curves, instead, 
correspond to a = 0.25, 0.5, 0.75, and 0.9 upon moving away from the diagonal line. 



(Note, however, that beyond the Gaussian approximation such a temperature depends 
upon the observable used to define it In equihbrium one has X°° = 1 (which 

is just a reformulation of the FDT) and P°° = fi. Instead, 7^ 1 is a signal of an 
asymptotic non-equilibrium dynamics and therefore we shall focus on this quantity for 
the dynamics we are presently interested in. 

Within the Gaussian approximation discussed in Sec. IS.ll the fiuctuation-dissipation 
ratio X can be easily calculated from the expressions in Eqs. ( 13 .Tp and ( 13. 8 p [see also 
Eqs. ( ]B.17p and ( lB.2ip ] for the response and correlation function, respectively, in terms 
of Mittag-LefHer functions: 

where we assumed t > t' and A = p^+r. For A ^ (e.g., far from the critical point r = 
or at criticality with p ^ 0) and long and well-separated times t,t',t — t' ^ (7/A)^/", 
one can easily see from Eq. ( 1B.19P that X~^ — )■ 1, confirming the expectation that the 
system equilibrates at long times, independently of the value of a > 0. On the other 
hand, for the fluctuation of the homogeneous mode p = at criticality one has A = 
and the FDR takes the simple form (originally derived in Ref. [28] for an anomalously 
diffusing particle) 

^^^io.c.it(^,0 = l+f^-lV " (4.44) 



t' 



which is a universal scaling function of the dimensionless scaling variable t'/t, reported 
in Fig. IH] for various values of a. In contrast to the white noise case a = 1, in the 
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presence of colored noise < a < 1, Xz}^ ^^.^^(t,t') does actually depend on t'/t and 
it interpolates continuously between the quasi-equilibrium regime t' ~ t, within which 
Xp-=o,crit — 1, and the non-equilibrium regime of well separated times t' <^ t, for which 
Xp-=o,crit — as it is generically observed in the case of coarsening dynamics [331 EH SSj • 
Beyond the Gaussian approximation, we can deduce an expression of the two- 
time dependent FDR and its limiting values from the scaling forms in Eq. fl4.22p 
and Eq. f l4.23p . First of all we note that for t ^ t' and p = one has dfC ~ 
^i2-r,)/z+e-a(^^ _ 0)t'^-'^~^Fci6,O) and R ^ t(2-''^^)/^(t/t')^i^i?(0, 0). We thus obtain 



X'- = hm hm (-] . (4.45) 

In the case d = 1 of dominant white noise this expression renders the well-known 
result X- = Ffi(0,0)/[(1 - e)Fc{0,0)]\a=i = X- [321 [H], i.e., X^ = 1/2 within the 
Gaussian approximation [^ W7\ HH] . The contribution of non-Gaussian fluctuations 
for D < 4 and up to the second order in the e-expansion have been calculated in [36] . 
in rather good agreement with Monte Carlo simulation (see Ref. [H] for a summary). 
Instead, if the colored noise dominates a = a < 1 and therefore the long-time limit 
X°° of the FDR in Eq. f l4.45p vanishes, formally corresponding to an infinite effective 
temperature as observed in coarsening processes. Note that this result holds at all orders 
in perturbation theory. Therefore, 

( X^ for a > aJD,N), 

[ for a < aciD,N), 

where both values do not depend on the actual value of a and therefore exhibits 
a discontinuity as a function of a upon crossing the line a = ac{D,N). Within the 
Gaussian approximation one can easily check the general result Eq. f l4.46p for X°°, on 
the basis of Eqs. f l4.2p and (14. 3p . Indeed the behavior of the correlation and response 
functions can be determined by taking limA_s.o ACo(0; A, k) and limx^o XRo(0] X, k), 
respectively, for the propagators at zero momentum and at criticality. It is then 
straightfoward to obtain X^ = 1/2 within region W and X^ = limx^o^K/^\ = 
within region C, which confirms our general results. Apparently, this result for X^ 
contradicts the corresponding one X^ = 1 for a freely diffusing particle in a super- 
Ohmic bath (corresponding to a > 1) found in [28], which our model reduces to within 
the Gaussian approximation. However, within the field-theoretical approach discussed 
here, it turns out that a super-Ohmic bath, responsible for a noise F with a > 1 in 
Eq. (12. 6p . is eventually controlled by the white-noise vertex and it is therefore unstable 
with respect to the effects of the interaction, which effectively generates such a vertex 
even though it was not present in the original coupling to the bath. Therefore, the 
white-noise result X^ = 1/2 does not only apply to the cross-over line ac{D,N) but 
it is valid within the whole region W. On the same footing, the results discussed here 
suggest that, at least in higher spatial dimensions, adding interactions to a system 
which displays super diffusion (corresponding to 2; < 2) results quite generically in a 
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sub-diffusive behavior {z > 2) as expected in the case of a diffusing particle (Gaussian 
approximation) with interactions. 

5. Summary and outlook 

We studied the purely dissipative critical dynamics of a model with an A^-component 
order parameter in D spatial dimensions, coupled to an equilibrium thermal bath which 
provides a colored thermal noise. We argued that the upper critical dimensionality of 
the model is = 4 and we used the framework of the field-theoretical e-expansion to 
account for the effects of non-Gaussian fluctuations in 4 — e spatial dimensions. 

Within the Gaussian approximation — valid for D > Dc — the equilibrium dynamic 
exponent z which controls the different scaling of space and time takes the values 



.(col) 



2/a for a < 1 



^0 = <^ " (5.1) 
^ = 2 for a > 1, 

where a characterises the algebraic long-time decay of the two-time correlation function 
of the noise, see Eq. (12. 6p . For a = 1 one recovers the white- noise result zj^^ = 2. The 
non-equilibrium 'initial slip exponent' 6, instead, vanishes. Depending upon the value 
of a the asymptotic long-time dynamics is effectively equivalent to one driven by white 
noise (Ohmic bath) for a > ac, wheras the effect of the colored noise is relevant for 
a < Oc- Within the Gaussian approximation o^c = 1, as demonstrated by the change in 
behavior of Zq given in Eq. (15. ip . 

In dimensions D < 4 the critical behavior is modified due to the relevance of the 
interaction term and of the non-Gaussian fluctuations. The value ac which controls 
the cross-over between the white-noise and the colored-noise dominated behaviors is 
modified by A^-dependent corrections of order and it therefore separates the two 
corresponding regions in the parameter space {a,D,N), named W and C in Fig. HJ 
respectively. The dynamical critical exponent z is given by 



^(coi) ^2^2 
a a 



, 4 1 

3 n 

3 2 



+ C(e^) within region C, 

(5.2) 

(AT + 8)2 r 3 2 J ^ ^^^'^ ^' 

The iV-dependent curve (13.32^ which separates regions W and C in the {a, -D)-plane is 
illustrated in Fig. IHfor iV = 1, 4, oo. Some comments are in order: 

(i) Upon decreasing D, the region W within which the Ohmic result is recovered 
extends beyond the Gaussian value ac = 1. 

(ii) The correction to the Gaussian value zo is positive within region W (2:0 = 2) and 
negative within region C (2:0 = 2/a). 

(iii) The exponent 2; is a continuous function of e and a: At the transition line between 
regions W and C one has z^^^ = z^^"^\ as can be easily verified by using Eq. (I3.32p . 
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(iv) In the large- TV limit the correction vanishes and the dynamic exponent z and Oc 
take their Gaussian values zq and = 1, respectively. 

For random initial conditions, i.e., with vanishing correlations and average order 
parameter, we determined the general scaling forms of the dynamic correlation functions. 
Within region C, such scaling forms differ from the ones valid in the presence of white 
noise only, studied in Ref. [10] and recovered within region W. We determined the 
corresponding initial-slip exponent 6 up to order 0{e) in the presence of colored noise. 
It is given by 



and the plot of the ratio between this value 6 and the reference 9a=i for the white noise 
is reported in Fig. [71 

In non-equilibrium conditions we also calculated the long-time limit X°° of the FDR 
for general a and N. The value of X°° in the presence of white noise is known analytically 
up to 0{e^) [36] and numerically via Monte Carlo simulations in various dimensions for 
models belonging to the universality class of the 0{N) model with dissipative dynamics 
(see, e.g., [Hj for a review). We proved that this result is recovered within region 
W. Instead, if the colored noise is dominant [a < ac{D,N)], i.e., within region C, 
we showed that X°° = 0. Therefore, the associated effective temperature is infinite, 
analogously to what is found in sub-critical coarsening [33l [34] . Our result for 
within the Gaussian approximation is only in partial agreement with the corresponding 
one derived in [28] for an anomalously diffusing particle — i.e., of a fractional Brownian 
motion — which our model reduces to within such an approximation. Indeed, in the 
presence of a super-Ohmic noise a > Oc = 1, one finds X°° = 1 [28] and super-diffusion 
z < 2 for the fractional Brownian motion, while we argue that X°° = X^ = 1/2 and 
normal diffusion z = z^'^ = 2 in our field theoretical model. This is due to the fact that 
even in the absence of a white-noise effective vertex in the original model, non-Gaussian 
fluctuations (induced by the interactions) generate it and turn it into the dominant one 
for a > ttc < 1 such that the white-noise result is recovered. 

In conclusion, noises correlated in time may affect significantly the equilibrium and 
non-equilibrium dynamical properties of systems close to critical points. In this respect 
it is important to note that the distinction between super-Ohmic {a > 1) and sub-Ohmic 
(a < 1) thermal baths does not fully correspond to having irrelevant (white) and relevant 
(colored) long-time correlations of the noise, respectively. Indeed, as shown in Fig. [U 
even a weakly sub-Ohmic noise with ac(-D, A^) < a < 1 is actually equivalent (in the 
RG sense) to an Ohmic (white) noise in the physical dimensions D = 3 and D = 2 
as far as the dynamical properties in the long-time limit are concerned. In addition, 
in the presence of interactions, a super-Ohmic bath does not result in a super- diffusive 
behavior {z < 2) but rather in the anomalous diffusion induced by the equivalent white 
noise, in contrast to what happens for the free fractional Brownian motion. 

The field-theoretical predictions for the relaxational Markov critical dynamics of 
systems belonging to the universality class considered here have been put to the 



6 



(iV + 2) 

4:{N + 8) 



d{a)TE{a)e + 0{e^) 



(5.3) 
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numerical test both via Monte Carlo simulations and by solving the Langevin equations 
with a variety of different methods (see, e.g., [15] and references therein). An instance 
of non-Markovian dynamics of the 0^-theory with a noise exponentially correlated in 
time was investigated in [19]. However, in this case one does not expect the long- 
time dynamics of the system to be affected by the finite memory of the noise. Dealing 
numerically with power-law correlated Gaussian noise is a significantly harder problem 
which remains basically open due to the difficulties in generating such kind of random 
process, see, e.g., [271 EO] ^ind references therein. 

One of the virtues of the approach we have followed here is that it can be easily 
applied to quantum critical dynamics [51]. For instance, the thermal bath can be 
modeled by a set of (quantum) harmonic oscillators coupled to all degrees of freedom 
of the system. Within the Schwinger-Keldysh formalism it is possible to derive a path- 
integral representation of the non-equilibrium dynamics [T8| [52] . Integrating out the 
oscillator variables one obtains an action similar to the one considered here [TH |52l |53] . 
The main difference with the classical case is that even Ohmic dissipation leads to 
retarded interactions. The present work is expected to provide at least a partial and 
preliminary insight into the more difficult problem of the analysis of quantum critical 
equilibrium and non-equilibrium dissipative dynamics |51] . 

Among other possible extensions of the present work, we mention the problem 
of understanding the effects of colored noise on sub-critical coarsening. The dynamic 
scaling hypothesis states that the late-stage phase ordering kinetics is governed by a 
length scale L{t) that, in models with no quenched disorder, typically grows in time 
as a power-law L{t) ~ A(T)t^/^'*. The dynamic exponent Zd (generically different for 
the dynamic exponent z at criticality) depends upon the kind of order parameter and 
the conservation laws [51] while the prefactor A(T) typically depends only weakly upon 
temperature T, is non- universal, and it vanishes upon approaching a critical point. (The 
matching with the critical growth is explained in [55].) In presence of colored noise this 
growth law might be modified, even though one usually expects thermal fluctuations 
not to affect the domain growth [5l]. We shall address this issue in a future study. 
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A. Fourier and Laplace conventions 

Within the present study we define the Fourier transform and its inverse via 

dte"*'"* and F{t) = — e'^' F{uj). (A.l) 

■oo J -oa 

Instead, for every A > the Laplace transform is defined as 

/■oo 

Fx= dte-^* F{t). (A.2) 
Jo 

In the main text we shall drop the hats, whenever this does not generate confusion. 



B. The equilibrium propagators 
B.l. Scaling in real time 

For a = 1 (white noise) the equilibrium propagators have a simple analytic form in 
the time domain [3l [TOl |35]. They can be calculated by applying an inverse Fourier 
transform to Eqs. (13 ■4p and (13.51) : 

i?o(p,t) = e(t)e-(p'+^)*/^- (B.l) 

Co(p,t) = ^r^e-(*''+")l*l/^". (B.2) 
+ r 

For general a dimensional analysis suggests that the critical (r = 0) Gaussian 
correlation Co with 7w = should scale as 

Co(p,t)=P"'/co(p'W77). (B.3) 
Using Eqs. (13. 3p and (13. 5p the equal-time correlator is given by 



Co(p,t = 0) = j ^C'o(p,c^) 



da; 
1 



27 sin(7ra/2) a; 


\a-l 




\uj\ 


|2" + 27(p2 + r)|a;| 


|°cos(7rQ;/2) + (p2 + r)2 



p'^ + r 



(B.4) 



Hence, we infer that /co(0) = 1- Naturally, we have /co(oo) = since correlations have 
to vanish in the long-time limit. Applying a Fourier transform to Eq. (IB. 31) it is easy to 
show that at criticality (r = 0) 

Co{x,t) = j^gc,hx'/\tn, (B.5) 

where the function gc^ reaches the asymptotic value 5^00(00) = r£;(D/2 — l)/(47r^/^). In 
order to deduce the leading behavior for gco{y) when u — )■ we start from the explicit 
expression of the noise kernel Ti^^ given in Eq. (13. 3p . After some algebra we obtain 

/■ d^p dw 2usin(7ra/2)|a;|"-ie*^+*?^-^ , , 

9cM= / 77^7^ ..21. .12^ , o..^2i..i.v„„.^„^,/oN , ^4 ' (B-6) 
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where u = 7x^/t" and z is an arbitrary unit vector. For a < 1 we neglect the 
contributions of in the denominator and we obtain 

' y 27r ' ' 7 (27r)^ (p2 + u|a;|"cos(7ra/2))2 ^ ' 

The integral over p is of In cos 7ra/2]) for D = 4 and the resulting integral 

converges for a < 1; consequently, 

gcoiu ^ 0) ^ 0{u In u). (B.8) 
By using FDT we derive 

Ro{x,t) = ^^S^. g'c,{ixyt'^)e{t). (B.9) 
In the white-noise case, the scaling function has the simple form 



r^p/2-1) 

= 4vr^/2 



1 _ + o(e) 

r,(f-i) 



(B.IO) 



with Te{s, x) = dy ^e *, whence we deduce for a = 1 and D = 4 

gc,{u->0) = O{u). (B.ll) 

For generic 7 and 7^ the scaling function is no longer a function of one variable. 
It is easy to show that 

Coix,t) = gcM^\V\tr,7MV\t\)- (B.12) 

Moreover, by using a similar argument as above one has for D = 4 

lim gco{u,v) = 0{{u + v)\nu), (B.13) 

t—>oo 

where u = 7|xp/t" and v = 7w|xp/t vanish with u/v finite. In the opposite short-time 
limit in which u and v diverge with u/v finite, 

limgc,{u,v) = Te{D/2 - l)/(47r^/2) (B.14) 

as for the purely colored problem. 

The equilibrium propagators can be written in terms of the generalized Mittag- 
LefHer functions Ea,b{z), as discussed in App. IB.2I 

B.2. Generalized Mittag-Leffter functions 
The Laplace transform of Ro{p, t) is given by 

where we defined A = + r and, in the case of colored noise, = 7A"^^. We formally 
expand this expression for small A: 

MP, A) = At , = ^ E (B-16) 

^ ^ 7A" 1 + A(7A°)-i 7A° ^ A"'^ ^ ^ 
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where the terms of the form l/A'^ (with Re/5 > 0) are recognized as the Laplace 
transform of Q{t)t^~^/rE{P), so that flB.16P is identified as the Laplace transform of 

= e(t)^E„,„(-Ar/7), (B.17) 

7 



where we have introduced the generalized Mittag-Lefiler function 

oo 

E^^p{z) = y^-—^—— with «,/3,2;eC, Re{a,/3}>0. (B.18) 

Note that this function reduces to an exponential for a = = 1: £'1^1(2) = e^, whereas 
for 2 G M Ha], 

- -) - - E f^pb^? + (B.19) 

The corresponding expression for the equilibrium Gaussian correlation function can 
be obtained from the FDT fl2.12p . Indeed, after integration Eq. fl2.12p takes the form 

r\t\ 

Co{p, t) = Co{p, t = 0) - / ds Ro{p, s) (B.20) 

Jo 

where we used the fact that, in equilibrium, C{x,t) = C{x,—t). Taking into account 
(1B.4P and the first line of flB.17p one readily finds 

Co{p,t) = ^E^{-A\trh) (B.21) 

where Ea{z) = Ea,i{z) is the Mittag-Lefiler function. 

The correlation function Co{p,t) in Eq. flB.2ip can also be expressed as the inverse 



Fourier transform of Co{p,u) reported in Eq. (13. 5p (see also Eq. (13. 3p ). After some 
suitable change of variables one finds the following scaling form 

Co{p,t) = jfcM\trh), (B.22) 

where 

2 r...„„.. v^-'sm{na/2) 



2 /"°° 

TT Jo 



y2a _|_ cos(7ra/2) + 1 



sin(7ra/2) f^^ _ cos(nV"t;V°) ^3) 



na/2 Jq _|_ 2t; cos(7ra/2) + 1 

is the explicit expression for the scaling function introduced in Eq. ( lB.3p . 



C. Calculation of S^'^ and £^'^ 

Starting from Eq. (I3.2ip we have for generic 7 and 7^ 



,2co,2, ,_,g'An{N + 2) 
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dt cos(at) £^-2^ gl^ 



t° ' t 



(C.l) 



The result of the integral in the first term in curly brackets is an analytic function of a 
that admits a Taylor expansion in powers of a^, i.e., 

Co + C2(J^ + Ci(J^ + • • • (C.2) 

with coefficients that, in principle, depend separately on 7, 7^ and I. The integral in 
the second term in curly brackets yields, instead, a non-analytic function of a that we 
can still express as a series: 

do + rfsc^' + . . . + rfsa-ifx^""' + . . . (C.3) 

where the term oc o"^""^ is due to the leading behavior of (7^^ for t — t- +00 [see Eq. ( IB. 81) ] 
which has to be subtracted for a < 1/3 in order to make the integral convergent at 
large t. If 3a — 1 > the limit a — can be safely taken and the white-noise vertex is 
renormalized by Cq -|- lio- If 5 011 the contrary, 3a — 1 < the contribution proportional 
to cr^"-^ is anyhow negligible (for a > 0) with respect to the term 70"" ^ which is 
already present in the tree-level vertex. Therefore, there is no renormalization of the 
colored- noise vertex and we can focus on the limit 7^ ^ 7, i.e., on the correction to the 
white-noise vertex only. Since we calculate evolution equations up to order we simply 
need to evaluate (]C.1|) m. D = A. We obtain 

^ g^AD{N + 2) 



n^f'^'^(0;7,7w) = 
X < zt-^ 
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dx X-'' gl^ 



7X^ 7„x^ 



+ 



dtr^ gl^ 



1^ Twf 
t° ' t 



(C.4) 



We are interested in the a ^ ac limit in which 7„ — )■ 00 and 2 = 2 -|- 0{e^\ By first 
using X I— >■ xij and t 1— )■ 7„£^/x^ we transform the two- variable scaling function into 
the one-variable white-noise one. Using then Eq. ( IB. 101) and ^4 = 27r^ we obtain the 
second and third line below. 



n2f0'2(0;7,7w) 



27„(72^z)(iV + 2) 



dxx ^g%^{^,x') 



+ 



/7w 



2-i^u^[N + 2) 
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/ dx 


"1 - e-"/^" 


/o 




4 





dxx 5-00(0'^ 



/x^ 



12 



In-. 

3 



(C.5) 
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Therefore, at the critical point, using the Wilson-Fisher fixed point value u* = 
6e/{N + 8) + C(e2) [37|, we find 

u^£0,2^0;^^^^)^u*'^^£'J = ^^^^^\n^e^ + O{e^ (C.6) 

We now compute S^'^ in D = 4. We start from Eqs. fl3.39p and fl3.40p . Using 
Eqs. (ID.SP and HB.lAh in the limit i ^ we obtain 

2ciUn ^ g'^MN + 2)71 -d [^dx 1 



144 X (27r)6 

u^{N + 2) 



(C.7) 



72 

Note that the term coming from the differentiation of Co{x, i^) in Eq. (13.391) with respect 
to ln£ vanishes in the limit £ — > [use Eq. fIB.Sp ]. At the critical point we obtain 



^ ^ 72 -2(iV + 8)2^ ^^-^^ 

D. Non-equilibrium propagators 

For a = 1 the Gaussian non-equilibrium propagators read in the momentum and time 
domain 

+ r I J 
Ro{p; t, s) = Q{t - s)e-(f'+'')(*-^)/^-. (D.2) 

For generic a, however, analogously compact expressions are not available and our 
analysis proceeds using Laplace transforms. In order to determine the response function 
Rq — consistently with the Gaussian approximation — we start with the linearized 
version of the Langevin equation Eq. fl2.2l) in the presence of an external perturbation 
h: 

"t 

dt' T{t - t')dt4{x, t') + (r - V^)0(f, t') = C(x, t) + h{x, t) (D.3) 



Calculating the expectation value of both sides with respect to the distribution of the 
noise eliminates the vanishing average {Q. The Laplace transform yields 

(AFa + / + r){$xmh = Hp) (D-4) 

in momentum space where we used the Dirichlet boundary condition t = 0) = 
[see discussion at the beginning of Sec. 14.1] . Note that the expectation value of the order 
parameter depends on h. The response propagator in the Laplace domain is given by 

^ 5.,. (D.5) 



Ro{p;X,K)6i 



{X + K){XTx+p^ + r) 
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The last equality follows from the fact that Shi{t)/Shj{s) = SijS{t — s) as a function of 
time translates into Sij/{X + k) in Laplace space, given that Jg°° dtdse~^^~'^'^5{t — s) = 

In order to deduce the correlation propagator we start directly from Eq. (ID. 31) with 
h = and we consider its Laplace transform: 

Up) = ,p (D-6) 

AFa + + r 

which yields 

{Ci,xCj,K) 



(ALa + p2 + r)(Kr« + p2 + r) 

5.r (D.7) 



(A + k) (AFa + p2 + r) (kF« + + r) 
In the last line we used the fact that dtdse-^*-'"^F(t - s) = (Fa + F«,)/(A + k). The 
propagators verify an 'initial time FDT'. We see that for a < 1 lim^^oo i^RoiP] -^j i^) = 
Ro{p, A) and lim«,^oo i^'^Co{p; A, k) = limfc_^oo '^^'""'^xRoip, A), with Ro{p, A) = 1/ (AFa + 
p"^ + r). In the time domain, the second identity reads 

dt'Co{p;t,t' ^0) r^t""^^ [ dsT{t- s)Ro{p;s,t' ^0). (D.8) 

Jo 

To derive this equation we used the convolution theorem for the Laplace transform C, 
that is C /g dt' f{t — t')g{t') (A) = C[f]{X)C[g]{X). In order to deduce the scaling of 
Eq. fID.SI) with respect to t' one observes that if XC[f{t)]{X) ~ A" for A — )■ oo then 
f{t) ~ t-" for t 0. 
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